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Abstract 

We propose an analysis of the time-optimal control of a dissipative two-level quantum system 
whose dynamics is governed by the Lindblad equation. This simple system allows one to use tools 
of geometric control theory and to construct its optimal synthesis i.e. to determine the set of all 
the optimal trajectories starting from a given initial point. We study different processes such as 
conversion of a pure state into a mixed-state and purification of a mixed-state. In particular cases, 
we show that dissipation is not undesirable and can help accelerating the control. 
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I. INTRODUCTION 



Manipulating quantum system by using time-dependent electric field remains a goal of 
primary interests in different physicalprocesses extending from the control of chemical reac- 
tions to quantum computing |4J . In recent years, active research has been performed 
to take into account the interaction of the system with the environment which represents 
more realistic situations but also more challenging control scenarios than for closed quan- 
tum systems 

aaa. 

Different control strategies have been proposed. In this paper we 
will not consider control techniques using non-unitary control such as measurement (see 



for instance the quantum Zeno effect) 



Q, Q, Q, Q 



or other strategies such 



as bang-bang pulses, strong coupling with another system or control that can actively act 
on the dissipation (see for a recent review and references therein). It has been shown 
that such methods are particularly efficient and can even halt decoherence. We restrict 
our study to unitary control. In this case and for the interaction with a markovian or a 
non-markovian bath, the fact that the control field cannot fully compensate the effect of 
dissipation largely enhances the difficulty of the control. This point has been rigourously 



shown in Refs. 



16 



17] for a dynamics governed by the Lindblad equation [18|, ]l9| . In this 



context, several studies using numerica 



control can still be achieved 120 



21 



22 



optimization techniques have proved that efficient 



23 



24 



25 



261 ] . Due to the complexity of realistic 



systems with multiple degrees of freedom, this purely numerical approach seems to be the 
only possible way to achieve control. Mo re g eometrical aspects of the control using mainly 



tools of geometric control theory 



27 



28 



291 ] can be formulated only for simplest quantum 



systems having few levels (typically two or three) or consisting in the coupling of spin 1/2 



particles. This has been done recently by a 
with closed 
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32 



33|, 13J, 135J, 136J, |37J, 



38 



arg e number of mathematical papers dealing 
39 ] or dissipative 4o| quantum systems. One 



of the objectives of the control has been the minimization of the total time of the process 
either with constraints 



30 



33] or no constraints on the laser intensity 34j. 



In this paper, we propose to do a step towards the geometrical analysis of the control 
of dissipative quantum systems by beginning with the simplest system possible, a two-level 



system governed by the Lindblad equation. We determine control fields which minimize 



the 



total time of the process by applying the Pontryagin maximum principle (PMP) 



27 



29 



41| . The maximum of the laser intensity is fixed to an arbitrary value. An increase 



2 



of this value leads to a reduction of the duration of the control. Note also that this cost 
functional seems particularly relevant in the context of a dissipative environment especially 
when the effect of dissipation is undesirable to reach the objective of the control. A second 
argument explaining our choice of cost functional is of mathematical nature. The resolution 
of the PMP is particularly simple in this case and can be done analytically. The geometrical 



description of the time-optimal control is also well developed especially on R 2 (see 28| for a 
recent overview) for affine systems with a drift term corresponding here to the dissipation. 
The main tools of this construction will be recalled throughout the paper. Finally, all this 
mathematical arsenal allows us to answer some physical questions such as the benefit that 
can be gained from dissipation. 

The paper is organized as follows. We first introduce the model for the Lindblad equa- 
tion of a two-level system. Writing the density matrix in the coherence vector form |42| and 
assuming some constraints on the control term, we restrict the problem to a control on a 
closed submanifold of R 2 with a true drift term i.e. an uncontrollable term which cannot 
be eliminated by a unitary transformation. This drift term is due to the dissipation. We 
consider different cases with different initial and target states that represent several situa- 
tions of physical interest. We formulate in Sec. II III the PMP with a time minimum cost 
functional. In Sec. \IV\ we analyze the structure of the reachable set from each initial state 
and we point out the particular role of the fixed point of the free-dynamics in this structure. 
Section |V] deals with the construction of the optimal syntheses by solving the PMP. We 
recall that the PMP allows us to derive extremal controls. Optimal trajectories are a subset 
of this set which can be determined in a second step by a direct comparison of different 
extremal trajectories or by geometrical arguments. Conclusions and prospective views are 
given in Sec JVIl Some technical calculations are reported in appendices \K\ and [B] 



II. THE MODEL SYSTEM 



We consider a dissipative two-level quantum system whose dynamics is governed by the 
Lindblad equation. The system is described by a density operator p(t) which is a positive 
semi-definite hermitian operator acting on a two-dimensional Hilbert space 7i. The evolution 
equation can be written as 

= [H + uH 1 ,p]+ijC D (p) , (1) 
3 



where Hq is the field-free Hamiltonian of the system, H% represents the interaction with the 
control field and Co the dissipative part of the equation. In the Lindblad equation, Cd can 
be written in the general form 

3 

(2) 
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where the operators L k are given by 
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7i2; 721 and T are positive real constants describing the interaction with the environment. 
7i2 and 721 correspond to population relaxations whereas T is the pure dephasing rate. We 
consider resonant fields, i.e. the frequency u of the laser is equal to the energy difference 
between the two levels. In the RWA approximation, the time evolution of p(t) satisfies the 
following Redfield form of the Lindblad equation 



(4) 



where u is the complex Rabi frequency of the laser field (the real and imaginary parts are 
the amplitudes of two orthogonal linearly polarized fields). V is the total dephasing rate 
which can be written as 
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(5) 



From Eqs. ([5]), we also notice that the requirement that Eq. (J4]) corresponds to a Lindblad 
equation implies the constraint T > ^(712 + 721) or equivalently T > 0. Equation (TJ]) is 
written in units such that h — 1. In the interaction representation, Eq. (j4]) becomes 

/ -Z7i2 -u* u «72i \ ( Pu\ 
-u —iY u P12 
u* — iY —u* p~2i 

\ «7l2 U* -U -2721 / \ P22 / 

Since Tr[p] = 1, the density matrix p depends on three real parameters which can be 

42| : Xl = 2M[p 12 }, x 2 = 2Q[p 12 ] and 
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given by the coordinates of the coherence vector 
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x 3 = P22 — Pii- From Eq. (jSj), one deduces that the coordinates Xi satisfy the following 
system of inhomogeneous linear differential equations 

x\ = -Txi + u 2 x 3 

x 2 = -Tx 2 - uix 3 , (7) 

^3 = (712 - 72l) - (712 + 72l)^3 + U!X 2 - U 2 Xi 

U\ and u 2 being two real functions such that u = U\ + iu 2 . As Tr[p 2 ] < 1, we also have 
x\+x^ + x 3 < 1 which defines the Bloch ball. The dynamics is called either unital if 712 = 721 
i.e. the fixed point of the free dynamics is the center of the Bloch ball or afline otherwise 



IS 



To simplify the study, we restrict the dynamics to a submanifold M C M 2 by assuming 
that the control field is real i.e. u 2 = 0. This hypothesis means that the control field is 
linearly polarized. With this choice of control, the last two equations of the system of Eqs. 
(J7J) are decoupled from the first one and the problem is reduced to IR 2 . The system of Eqs. 
(J7|) then becomes 

x 2 = -Tx 2 - ux 3 
x 3 = 7- - 7+^3 + ux 2 

where 7_ = 7 12 — 721 and 7+ = 712 + 721- The coordinate X\ is set initially to 0. To simplify 
the notation, the index 1 of U\ has been omitted when confusion is unlikely to occur. Note 



that the analysis of the optimal control on R 3 is considerably more complex 29( and goes 
beyond the scope of this paper. Equations (jHJ) can be written in a more compact form 

x = F + uG , (9) 

with the vector x of coordinates (x 2 , x 3 ) and the two vector fields F and G defined by 

(10) 



1 -tx 2 y 


and G = \ 






I 7- - 7+^3 J 







Finally, a straightforward calculation shows that the fixed-point of the free-dynamics is given 
for 7 + 7^ by 

' x 2 = 



7- 

X 3 = — 

7+ 



(11) 



which is therefore either afline or unital according to the value of 7_. The irreversibility 
of the dissipation effects is reflected in the fact that the vector field F is a true drift term 



which cannot be eliminated by unitary transformations 
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III. PONTRYAGIN MAXIMUM PRINCIPLE 



We analyze the optimal control of this two-level system with the constraint of minimizing 
the total time of the control. We assume that the field u is bounded by 

\u\ < 1 . (12) 

Equations ([8]) being linear, other bounds for u can be considered from a standard rescaling 
of the time and the dissipative constants T, 7+ and 7_. The Pontryagin maximum principle 



3, 



28 



29 



4l| is formulated from the following pseudo-Hamiltonian Hp 

H P = p-(F + uG)+ Po , (13) 



where p = (^2,^3) £ (M 2 )* is called the adjoint state and po is a negative constant. We recall 
that the cost function for the minimum time problem is equal to 1. This term is multiplied 
by the constant po i n the Hamiltonian Hp. The Pontryagin maximum principle states that 
the extremal trajectories maximize Hp i.e. 

H max (x, p) = Max\ u \<iH P (x, p,u) . (14) 

The coordinates of the extremal vector state x and of the corresponding adjoint state p 
fulfill the Hamilton equations 

x 2 = —Txo — UX-x 

(15) 

x 3 = 7_ - 7+0:3 + ux 2 
and 

p 2 = Tp 2 -up 3 

(16) 

p 3 = 7+p 3 + up 2 
where u is here the extremal control given by Eq. ( Fl4l . 

We can now pass to the construction of the optimal syntheses. The construction begins 
with the introduction of two sets of points A^ 1 (0) and A^ x (0) denoted Ca and Cp which 



divide M in different regions 28j]. and A# are two functions from M to R defined as 
follows 

' A A (x) = Det(F,G) 
A B (x) = Det(G,[F,G}) 



where Det is the determinant of two vector fields and [., .] their commutator. In our case, 
simple algebra leads to 

Aa(x) = -Vx\ + 7-^3 
A B (x) = 2Tx 2 x ?J - 2<y + x 2 x 3 + 7_x 2 

The major role of these two sets in the resolution of the optimal control problem will be 
detailed in Sec. M We can already say that Ca and Cb are responsible for qualitative mod- 



ifications of the optimal trajectories 28|, |29(. A preliminary step thus consists in analyzing 
the structure of the sets Ca and Cb when T, 7+ and 7_ vary with the conditions 7+ > 
and T From a formal point of view, this can be done by introducing the feed-back 



group [29|] but the dynamics being here bilinear [Eqs. (jHj)], we can directly compute the sets 
Ca and Cb and their relative positions. 

If T 7^ 7 + then the set Cb corresponds to the union of the two lines 

x 2 = , (19) 

and 

^3 = ~ 7 ~ , (20) 

with the restriction that \x^\ < 1. In the case r = 7+, Cb is only composed of the vertical 
line of equation x 2 = 0. The points (x2,x 3 ) of the Bloch ball that are solutions of the 
polynomial equation 7+Xg — 7_x 3 — Yx\ = belong to Ca- In the case 7_ 7^ 0, Ca is 
therefore the union of two parabolas. This set is either above or below the line £3 = 2r "l 7 2 T y+ 
according to the signs of 7_ and of T — 7+. For 7_ = 0, this set is reduced to the origin of 
the Bloch ball. Figure [T] displays the submanifold M together with the sets Ca and Cb for 
a given value of the parameters. We finally notice that the fixed-point of the free dynamics 
and the center of the Bloch ball belong to Ca- 

IV. REACHABLE SETS AND CONTROLLABILITY 

We consider four different qualitative cases of control which allow one to study several 
physically relevant situations : 

• Case (a) : Conversion of a pure state into a mixed state with a unital Lindbladian 

(7-/7+ = 0, r> 7+ + 2). 
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FIG. 1: Division of the manifold M by the sets C A = A^ x (0) (in solid line) and C B = A^O) (in 
dashed and dot-dashed lines) for 7_ ^ 0. The dashed and dot-dashed lines represent respectively 
anti-turnpike and turnpike arcs which are introduced in Sec. |VJ The exterior circle in solid line 
corresponds to the limit of the Bloch ball in the plane {x2, X3). The small open circle indicates the 
position of the fixed-point of the free-dynamics. Numerical values are taken to be V = 1, 712 = 0.1 
and 721 = 0.3. 

• Case (b) : Conversion of a pure state into a mixed state with a unital Lindbladian 

(7-/7+ = 0, 7+-2<r< 7+ + 2). 

• Case (c) : Purification of the completely random mixed state which corresponds to 
the center of the Bloch ball (7-/7+ = — 1, L > 7+ + 2). 

• Case (d) : Conversion of a pure state into a mixed state with an affine Lindbladian 

(7-/7+ = -0.5, r> 7+ + 2). 

The numerical values we have chosen for illustrations are given in Table H Although the 
TABLE I: Numerical values of the dissipative constants in arbitrary units. 
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choice of the parameters of Table d will become clearer in Sec. |Vj some comments can 
already be made. This choice both depends on the structure of the sets Ca and Cb and on 
the characteristics of two particular extremals denoted X— and Y — which start at the initial 
point, and correspond respectively to a constant control equal to -1 and 1. As detailed in the 
appendix [Bj the X— and Y — trajectories are either pseudo-periodic or aperiodic according 
to the sign of the discriminant A = (T — 7+) 2 — 4 of the system of Eqs. ([8]). An exact 
resolution of the dynamics shows that the trajectory of the system is aperiodic if A > and 
pseudo-periodic otherwise. This point is summarized by the diagram [10] of the appendix 
iBl In Table H] we have chosen for three of the four examples T such that T > 7+ + 2 to 
simplify the local structure of the optimal synthesis around the fixed-point of the dynamics. 
A pseudo-periodic trajectory is locally a spiral in the plane (x 2 ,x 3 ) around this fixed point 
which makes the analysis more complex (see Sec. |V|) . 



A. Purity and limits of the dynamics 

Before analyzing the reachable set and the controllability of each example, we begin by 
some general comments about the dynamics. 

We first show that the field cannot locally compensate the effect of dissipation. The 
purity of the quantum state is defined by the function 2Tr[p 2 ] — 1 = x\ + x\. Pure states 
are thus on the unit circle of the (X2, £3) plane. Simple algebra then leads to 

= -Txl + 7_* 3 - l + xl , (21) 
and we notice that this derivative does not depend on u, which completes the proof. A more 



general proof of this point is given in Refs. 16|, [l7| . We also point out that for points x 
where A^(x) < then dT ^ p - < and inversely if A^(x) > then dT ^ p - > 0. The curve Ca 
divides the plane (x2,Xs) into a region where the purity of the state locally increases and a 
region where it locally decreases. On the boundary Ca, the purity is preserved. This point 
can be qualitatively understood as follows. We recall that the purity is equal to the square 
of the distance to the origin. The conservative vector field G is orthoradial (i.e. normal to 
radial vectors) for each point (^2,^3) 7^ (0,0) of the manifold. The dissipative vector field 
F does not modify the purity of the state if the radial component of F vanishes i.e. if F is 
parallel to G which is the definition of the curve Ca- 



9 



We next analyze the fixed points of the dynamics when the field is on, which are defined 
by F + uG = 0. Since F and G are parallel, the fixed points belong to the curve Ca- The 
field- free limit point is the point of this line of maximum purity (see for instance Fig. [T]). 
This shows that the dissipation alone allows to reach the state of maximum purity. Inversely, 
one can ask if every point of the curve Ca corresponds to a limit point of the dynamics. The 
answer is positive for a real non-bounded control u since the limits can be written 



— «7_ 
7- 

7++M 2 / r 



(22) 



B. Reachable sets 



3 . 



We begin by recalling some results of Refs. [la, [l7| about controllability of dissipative 
systems. A quantum dissipative system of finite dimension whose dynamics is governed 
by the Lindblad equation is generically accessible but not controllable. The accessibility 
property characterizes the fact that the system can be driven in every direction of the state 
space. Moreover, the concept of accessibility does not take into account the reversibility or 
the irreversibility of the process. The lack of controllability is measured by the non small- 
time controllability of the system which illustrates the irreversibility of the dynamics. This 
kind of system is not small-time controllable because the field cannot locally compensate 



the effect of dissipation as shown by Eq. (IzTil . The accessibility property can be checked 
by the computation of the dimension of the dynamical Lie algebra L of the system 

HQ 

We introduce for that purpose the density matrix p of coordinates (1,22? ^3) which is given 
in the basis of the coherence vector and we rewrite Eqs. (JSJ in matrix form as follows 



p = C F p + uC G p 
where Lp and Cq are 3x3 matrices given by 



(23) 
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(24) 



L is the Lie algebra generated by Cf and Cq- If T 7^ 7+, a direct computation shows that 
L is either isomorphic to gt(2) or to the semi-direct sum g[(2)(s)IR 2 of respective dimensions 
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4 and 6 for 7_ = or 7_ ^ 0. The system is therefore accessible for r ^ 7 + . 

We now determine the reachable sets from their respective initial states of the four exam- 
ples. The reachable sets can be constructed by the explicit construction of all the trajectories. 
We denote by 7£(x ) the reachable set from x . We first search for in this section the bound- 
ary of the reachable sets. Then in Sec. |V]we show that all the points inside the boundary 
are attainable. We consider for that the X— and Y — trajectories starting from the initial 
point and from the field-free fixed point of the dynamics. A qualitative change occurs when 
the X— or Y — trajectories cross Ca since the angle between the vector F(x) and G(x) 
changes its sign. The determination of all the extremal trajectories in Sec. |V] allows us to 
complete the construction. 

We apply these remarks to the cases (a) and (b) where the initial point has coordinates 
(x 2 = 0,x 3 = 1) and Ca is reduced to the origin (x 2 = 0,x 3 = 0). The origin is attained 
asymptotically by the dynamics when t — > +oo. The Y — and X— trajectories, which have 
the same initial and final points, are therefore global boundaries of the reachable set. This 
point is illustrated in Figs. Eh and [2b for the cases (a) and (b), where the reachable set 7Z 
is in grey. 

For the case (c), the X— and Y— trajectories starting from (x 2 = 0,^3 = 0) intersect 
asymptotically Ca- We also consider the X— and Y — trajectories starting from the field-free 
fixed point which define two new regions. As an infinite time is necessary to reach this point, 
these two trajectories do not belong to the reachable set. 7Z is the union of all these regions 
as shown in Fig. [2b. We will check in Sec. |V]by constructing the extremal trajectories that 
every point of this set is effectively attainable. 

For the case (d) the situation is more difficult to analyze as the X— and Y — trajectories 
starting from (x 2 = 0,^3 = 1) intersect transversally Ca- From this point of intersection, 
the two trajectories do not correspond anymore to the boundary of 7Z. After the crossing of 
Ca, the X— or Y — trajectories originating from the preceding curves define two new regions 
of 1Z. We finally consider the limit point of the field-free dynamics which can be attained as 
t — > +oo and the X— and Y — trajectories starting from this point. These two lines define 
two new boundaries which intersect asymptotically Ca but which do not belong to TZ. TZ 
is therefore composed of all the regions constructed above and is displayed in Fig. [2H. As 
before, only the determination of the extremal trajectories will complete the construction. 
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FIG. 2: Reachable sets 1Z (in grey) for the cases (a), (b), (c) and (d). The sets Ca and Cb are 
represented in dashed lines. The position of the field-free fixed point is indicated by the small open 
circle. In (c) and (d), the trajectories starting from the field- free fixed point do not belong to 1Z 
(see text). 

V. TIME- OPTIMAL CONTROL 

The application of the Pontryagin maximum principle is particularly simple at least 
locally in the case of time minimum control. The goal of this section is to solve this problem 
by constructing the optimal synthesis i.e. the set of all the optimal solutions starting form 
a given initial point xo and reaching a point of TZ(xq). For time optimization, the optimal 
controls are composed of piecewise constant parts (u = ±1) and of singular controls. Some 
preliminary work in Sec. IV Al has to be done before determining the optimal solutions. 
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A. Preliminary 



We define two functions $ and 9 28]. $, called the switching function, is given by 



$(t) = p.G = -p 2 x 3 + p 3 x 2 . (25) 

Using Eq. (fT3|) . standard considerations of maximization of the pseudo-Hamiltonian Hp lead 
to the conclusion that the extremal field at time t is given by u = sign(Q(t)) if ^ 0. A 
time to such that the control changes sign (i.e. such that $(to) = 0) is a switching time. If 
$ vanishes on an interval [to^i] then the corresponding trajectory is called singular and is 
referred to as a Z— trajectory. An important point is the fact that the Z— trajectories lie 
in the set Cb and that the corresponding control u = <fi is singular i.e. different from 1 or 
-1. <p can be calculated by imposing that ^A B (x(t)) = on [io^i]- In our case, it can be 
shown that 

^ (x(t)) = = ^7 X2 + ^ 3 - (26) 

One arrives after simple algebra to 

-x 27 „(r-2 7+ ) -2x2X3^1- T 2 ' 
2(r- 7+ )(x|-x|)-7_x 3 



<M X ) = ^7F 7~T1~j1 3x — • ( 27 ) 



For the line x 2 = of Cb, this leads to 0(x) = 0. For the line £3 = 2 (r- 7+ ) 1 we obtain 

0(x) = . (28) 

2(r - 7+)x 2 

The control is admissible if |0(x)| < 1 which implies here the condition 

,0-^-2^), (29) 

21 - 1 2(r- 7+ ) 

Not every Z-trajectory can be an optimal trajectory. We nrtrodnce to characterize this 
point the notion of turnpike and anti-turnpike curves [28]. Let x ^ Ca U and the 
function /(x) = — ^(x) - ^ turnpike curve is an arc lying in Cb such that for every point 
x of this arc A^(x) 7^ and A(x) and Y(x) are not tangent to Cb- It is also assumed 
that A(x) and Y(x) point to opposite sides of Cb which define two regions fl x and Q y . If 
/(x) > (resp. /(x) < 0) on Q y and /(x) < (resp. /(x) > 0) on Q x then the arc is a 
turnpike (resp. anti-turnpike) arc. The relation with the optimal synthesis can be stated as 
follows. Using for instance the clock form a (see appendix |A|), it can be shown that the anti- 
turnpike trajectories are not optimal. Figure [I] displays the turnpike and the anti-turnpike 
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curves for particular values of the parameters. 

The switching function $ is a powerful tool that gives all the informations on the extremal 
trajectory. However, a global analysis of the extremals using only $ is difficult because all 
the initial values (p20;P3o) have to be tested. A more global point of view is given by the 
function 9 which has the advantage not to depend on p. 

Let v be the vector of coordinates (x 2 , £3). Deriving Eqs. (jHJ) with respect to time, it can 
be shown that v satisfies the following system of equations 

= — Yvo — uv 3 

(30) 

= -7+^3 + uv 2 

We next introduce the vector v{t). At time t, let us assume that the dynamics reaches the 
point (x2{t), Xs(t)) and that the vector field G at this point is given by G{x2(t), x%(t)) = G(t). 
v(t) is defined as the solution at time of Eqs. (13"01 that at time t is equal to G(t). The 
dynamics is thus propagated backwards during the time t. 9 is then defined as 

9(t)=arg(v(0),v(t)) , (31) 

where the angle is measured counterclockwise. The definition of 9 originates from the fact 
that the product p(t) ■ v(t) is a constant j^]. We thus have 

$(t)=p(t)-G(t)=p(0)-v(t) . (32) 

The function 9 allows to determine some properties of the extremal trajectories. We briefly 
recall these points here. The reader is referred to [281 for rigorous definitions and complete 
proofs. It can first be shown that 

sign{9) = sign^s) , (33) 

which means that the zeros of 9 are located onC^. A switch can occur if 9 > and 9 > or 
9 < and 9 < 0. In addition, the variation of 9 between two switches or between a switch 
and a singular control is equal to modulo it. This latter point can be easily understood 
from Eq. (132]) . Other properties of 9 will be detailed during the construction of the optimal 
synthesis. 
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B. Time-optimal syntheses 



The parts IV B 11 IV B 2\ IV B 31 and IV B 41 are rather technical and describe the way to 
obtain the optimal syntheses. Section IV B 51 details conclusions on the role of the dissipation 
that can be gained from the resolution of the time optimal control. 

1. Case (a) 

In the case (a), the optimal trajectories are either bang or bang-bang. A bang trajectory 
is a trajectory associated to a single value of the control u = +1 or u = —1. A bang-bang 
trajectory is the concatenation of an X— and an Y — trajectories. We denote by X * Y 
such a concatenation where the Y — trajectory comes first. Here the maximum number 
of switches is thus equal to 1. The elimination of extremal trajectories selected by the 
Pontryagin maximum principle has been done through the clock form a (see appendix |A] 
for details) and the symmetry of the diagram with respect to the line x<i = 0. We recall 
that the clock form can only be used for trajectories belonging to one of the four quadrants 
defined by Cb and which do not cross Ca- Figure [3] displays the optimal synthesis for this 
problem and the evolution of the angle 9 along the Y — trajectory starting from the initial 
point (0, 1). The plot of 9 shows that the switches are always permitted along the X— or 
the Y — trajectory. Note also that the form of the curve representing 9 is related to the sign 
of A. More precisely, if T > / y + + 2 then 9 is a monotonically decreasing function and if 
T < 7+ — 2 (case not treated here), 9 starts increasing, passes through a maximum when 
the trajectory crosses Cb and then decreases. The intermediate case 7+ — 2<r<7 + + 2 
corresponds to the case (b). Using the clock form, it can be shown that only one switch is 
possible for a trajectory in a given quadrant. This point can also be determined from more 



general considerations detailed in Ref. [281 ] . The rest of the optimal synthesis is deduced 
from the symmetry of the problem. The line 22 = is called an overlap curve denoted by 
K as it is the locus of points reached by two optimal trajectories. We finally notice that the 
dissipation alone is not used here by the control and therefore cannot help accelerating the 
control. We recall that the singular control on the vertical line oi Cb is given by u = 0. 
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FIG. 3: The top figure represents the optimal synthesis for the case (a) (see text). The dashed 
lines indicate the locus Cb and the small open circle the fixed point of the dynamics. The bottom 
figure displays the evolution of 9 as a function of the time for the Y— trajectory starting from the 
point (0, 1). 

2. Case (b) 

The situation is a little more complex in the case (b). For x% > 0, the synthesis is similar 
to the one of case (a) i.e. the trajectories are either bang or bang-bang. The function 9 is here 
periodic, the times where 9 vanishes correspond to the crossing of Cb- This function tells 
also us that the initial Y — and X— trajectories cannot switch in their respective second 
quadrant (i.e. the second quadrant they go through) since 9{t) > and 9{t) < 0. By 
symmetry with respect to the line X2 = 0, one deduces that theses curves are optimal up 
to this line. A singular control along the horizontal line of Cb is optimal from the point 
of intersection between the initial X— or Y — trajectories and this line. Since 7_ = 0, the 
singular control is given by <fi = 0. From this singular line originate optimal trajectories 
with control u = ±1. Using the clock form a, it can be shown that these trajectories can 
not switch again. To summarize, in this second part we have constructed optimal controls 
of the form Y * Z * X , Y * Z * Y , X * Z * X and X * Z * Y. This optimal synthesis is 
represented in Fig. HI Here, one sees that the dissipation alone (with u = 0) accelerates the 
control to reach a point near the origin only along the horizontal direction and for 23 = 0. 
We finally give in Fig. [5] an example of the comparison of two extremal trajectories. The 
time to reach respectively the points A and B from the initial point (0,1) by the Y — and 
X— trajectories is the same. We would like to attain the point C from A or B. We use 
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FIG. 4: Same as Fig. [3] but for the case (b). The dotted-dashed line represents the singular 
trajectory S. The small insert is a zoom of the optimal synthesis near the origin. 




FIG. 5: Use of symmetry to determine the optimal trajectory (see text). 

either the X— trajectory from B or a concatenation of a Y — and a Z— ones from A. The 
clock form a cannot be used here since the trajectories belong to two different quadrants. 
We therefore consider the symmetric image AD of the curve BD with respect to £2 = 0. 
The two extremals are thus in the same quadrant and a can be used. We conclude that the 
curve Y * Z * Y is optimal to reach the point C. 

3. Case (c) 

In this case, the singular line X2 = is optimal. The singular control is equal to zero. 
Switches can occur from the initial X— and Y — trajectories but they do not lead to optimal 
trajectories. Inversely, X— and Y — trajectories originating from S are found to be optimal. 
When two extremals cross Ca, a cannot be used and a direct numerical comparison is then 
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FIG. 6: Same as Fig. d]but for the case (c). 

performed. In this case, the function 9 tells us that these curves do not switch. Trajectories 
originating from S that have a switch are therefore not extremal. We finally point out that 
the dissipation is beneficial for the control and decreases the duration needed to purify the 
system. The optimal controls are here of the form X—, Y-, X * Z or Y * Z. Figure [6] 
displays the optimal synthesis for this problem. 



4- Case (d) 

The case (d) is the more complex one and corresponds roughly to the composition of 
cases (b) and (c). The difficulty lies in the global structure of the control or, in other words, 
in gluing the two preceding local analysis. For x 3 > 0, the optimal synthesis is similar to the 
cases (a) or (b) with bang or bang-bang trajectories. The bottom of the optimal synthesis 
from the point of intersection of the initial X— and Y — trajectories is similar to the case 
(c). 

We now describe the central part of the synthesis. The horizontal singular line of Cb does 
not correspond to a singular trajectory since |0(x)| > 1 which is a non- admissible control. 
Consider the first point of intersection of Cb with the X— or F— trajectory. Following 2s|, 
we know that a switch curve C originates from this point. To determine the exact locus of 
C, we use the 9 function since this function varies by modulo it between two switches. C 
is constructed numerically. By definition of C, every trajectory which crosses C switches on 
C. We have found that C, Cb and Ca intersect in one point, the origin O. 

Since the line x 2 = is turnpike for 2= < x 3 < 0, we can ask if this singular trajectory is 
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FIG. 7: Same as Fig. 0]but for the case (d). The dotted line represents the switch curve C. 



optimal i.e. if we can have a local optimal synthesis of the form given by Fig. [HJ To answer 
this question, we use the switching function $. For x(£) G C U S, $(£) = i.e. the vectors 
p(i) and G(pc(t)) are orthogonal. Since the direction of G is known (G is orthoradial), one 
can deduce the direction of p(£). Let z x and z 2 be two points belonging respectively to C 
and S. The vectors G(z) associated to these points are schematically represented in Fig. 
[HJ We let now the states z x and z 2 go to (0, 0) and we determine the directions of the 
different adjoint states. We recall that the Pontryagin maximum principle states that p 
is a continuous function which does not vanish. When z x goes to (0,0), one deduces by a 
continuity argument that p 1 is vertical in O. When z 2 goes to (0,0), the limit direction of 
p 2 is given by the switch curve C. To respect the continuity of p, one sees that C has to 
be tangent to the line x 2 = in O. Due to the complexity of analytical calculations, we 
have checked numerically that this is not the case. The singular line for x 3 < is therefore 
not optimal. The trajectories of the form Y*X*Y or X*Y*X are thus optimal up to 
x 2 = 0. In addition, when the initial X— and Y— trajectories cross the angle between 
the vectors F(x) and G(x) changes its sign. New optimal trajectories originate from this 
point of intersection and correspond to two new regions of the reachable set. 

5. Qualitative conclusions on the dynamics 

From the results obtained in the preceding sections, some qualitative conclusions can be 
made with respect to the dissipation effect on the time optimal control of the dynamics. 
The dissipation is not undesirable when the dissipation allows to purify the system [cases 
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FIG. 8: Possible optimal synthesis around the origin. 



(c) and (d)] and help accelerating the control. In contrast for the conversion of a pure state 
into a mixed state, the dissipation alone increases the duration of the control and its effect 
is not beneficial for the control [cases (a) and (b)]. The example (d) summarizes well the 
situation. As long as the purity of the state decreases, it is advantageous to use a control 
field but when the purity starts increasing the dissipation alone becomes more efficient. 

Finally, we point out that a visual inspection of the optimal syntheses allows to check 
the form of the reachable sets given in Sec. IIVBI 



VI. CONCLUSION 



We have investigated the time-optimal control of a two-level dissipative system. For 
different initial and final points, we have constructed the optimal syntheses which allow us 
to rigorously conclude on the beneficial role of the dissipation in order to reduce the duration 
of the process. An open question is the generalization of the present study to more complex 
quantum systems with for instance two control parameters or three levels. The determination 
of the optimal syntheses for these systems requires more sophisticated mathematical tools 
than those used in this paper [29 1. In addition, such an analysis could also be performed 
on other problems of quantum control. We point out that the question of measurement in 
the context of geometrical control theory is particularly interesting. Measurements can be 
viewed as a non-unitary control which allows to enlarge the reachable sets. For instance, 
we can consider systems not completely controllable by unitary control but which become 
controllable if measurements are used. A simple example is given by a finite number of 
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levels of the harmonic oscillator with a constant dipolar interaction only between adjacent 



energy levels 



43]. From a formal point of view, this problem raises the challenge of the 



precise formulation of the Pontryagin maximum principle in this context which is to our 
knowledge an open question. We notice that all these geometrical tools would allow us to 
answer some physical relevant questions such as the time at which measurements have to be 
done to minimize the duration of the control. 



APPENDIX A: THE CLOCK FORM 

We derive in this section the expression of the clock form denoted a j^]. By definition, 
the clock form is a 1-form which fulfills the following conditions 

= 1 • (Al) 
= 

A solution of this system exists except on the set Ca where F and G are collinear. If we 
write a as a = a 2 dx 2 + ot^dx^ then simple algebra shows that a 2 and 0-3 are solutions of the 
system 




a 2 (-Tx 2 ) + 013(7- - 7+^3) = 1 
a 2 x 3 = a 3 x 2 



(A2) 



We obtain that 



_ -2'2 

^1-7-^3+7+^3 

_ -S3 

r^|-7-^3+7+^3 

From the 1-form a, we can define the 2-form da which is given by 



(A3) 



, . da 3 da 2 s , Al , A . , 

da = (— — — )dx 2 A dx% , (A4J 

ox 2 (JX3 

and reads after some calculations 

2rx 2 x 3 + 7_x 2 - 2^ + x 2 x 3 

da = — r— 5 5T- — dx 2 A dx 3 . { A5) 

[Tx 2 2 - 7_x 3 + 

If we write da as da = g(x 2 , x 3 )dx 2 A dx$ then one sees that g(x 2 , x 3 ) = on Cb and that 
the function g has a constant sign in the regions delimited by the lines of Cb- This point is 
displayed in Fig. [91 As suggested by its name, the clock form is a form which allows one to 
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FIG. 9: Sign of the function g in the plane (a^,^)- The large solid lines indicate the position of 
the set Cb- 

determine the time taken to travel a path and to compare the extremals. Let 7 be a path 
in the plane {x 2 ,x^) which does not cross Ca and T the time of travel along 7. We have 



ft(x)<it 



a(F)dt = T 



(A6) 



'0 jo 

We consider now two paths 71 and 72 starting and ending at the same points and respectively 
associated to the durations T\ and T 2 . One shows by using the Stokes theorem that 



7i-T 2 



a: 



a: 



da 



(A7) 



'71 ^ 72 J D 

where D is the surface delimited by 71 U —72- For paths 71 and 72 which lie in one of the four 
quadrants defined by Cb, we can straightforwardly determine the time-optimal trajectory. 



APPENDIX B: ANALYTICAL DETERMINATION OF THE DYNAMICS 



In this section, we solve analytically the dynamics of the system given by the following 
system of differential equations 



(Bl) 



x 2 = —Tx 2 — ux 3 
X3 = 7- - 7+^3 + ux 2 

We assume that u = ±1 and we denote by u = e the control term. Combining Eqs. (IB II) . 
one arrives to equations that depend on only one variable x 2 or x 3 



x 2 + (r + 7 + )a? 2 + (1 + r 7+ )x 2 + e 7 _ = 

fa + (r + 7+)^ 3 + (1 + r 7+ )x 3 - r 7 _ = 



(B2) 
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FIG. 10: Sign of the discriminant A as a function of the parameters T and 7+. The parameters 
are in arbitrary units. The zone in grey is excluded for Lindblad dynamics. 

Equations (1B2j) are linear inhomogeneous differential equations of second order whose dis- 
criminant A is given by 



A qualitative change of the solutions is expected according to the sign of A, i.e. a transition 
from a pseudo-periodic solution for A < to an aperiodic solution for A > 0. This point is 
summarized in Fig. [TUJ Note that the parameter 7_ plays no role in the structure of these 
trajectories but modifies their limits points. Simple algebra allows one to calculate the exact 
solutions. For instance for A > 0, we obtain 
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Time optimal control of a two-level dissipative quantum system 
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We propose an analysis of the time-optimal control of a dissipative two-level quantum system 
whose dynamics is governed by the Lindblad equation. This simple system allows one to use tools 
of geometric control theory and to construct its optimal synthesis i.e. to determine the set of all 
the optimal trajectories starting from a given initial point. We study different processes such as 
conversion of a pure state into a mixed-state and purification of a mixed-state. In particular cases, 
we show that dissipation is not undesirable and can help accelerating the control. 

PACS numbers: 32.80.Qk,03.65.Yz,78.20.Bh 



I. INTRODUCTION 

Manipulating quantum system by using time- 
dependent electric field remains a goal of primary inter- 
ests in different physical processes extending from the 
control of chemical reactions [? ? ? ] to quantum com- 
puting [? ]. In recent years, active research has been per- 
formed to take into account the interaction of the system 
with the environment which represents more realistic sit- 
uations but also more challenging control scenarios than 
for closed quantum systems [???]. Different con- 
trol strategies have been proposed. In this paper we will 
not consider control techniques using non-unitary control 
such as measurement (see for instance the quantum Zeno 
effect) [???????], or other strategies such as 
bang-bang pulses, strong coupling with another system 
or control that can actively act on the dissipation (see 
[? ] for a recent review and references therein). It has 
been shown that such methods are particularly efficient 
and can even halt decoherence. We restrict our study to 
unitary control. In this case and for the interaction with 
a markovian or a non-markovian bath, the fact that the 
control field cannot fully compensate the effect of dissi- 
pation largely enhances the difficulty of the control. This 
point has been rigourously shown in Refs. [? ? ] for a 
dynamics governed by the Lindblad equation [? ? ]. In 
this context, several studies using numerical optimization 
techniques have proved that efficient control can still be 
achieved [???????]. Due to the complexity of 
realistic systems with multiple degrees of freedom, this 
purely numerical approach seems to be the only possible 
way to achieve control. More geometrical aspects of the 
control using mainly tools of geometric control theory [? 
? ? ] can be formulated only for simplest quantum 
systems having few levels (typically two or three) or con- 
sisting in the coupling of spin 1/2 particles. This has 
been done recently by a large number of mathematical 
papers dealing with closed [?????????? 
] or dissipative [? ] quantum systems. One of the ob- 
jectives of the control has been the minimization of the 
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total time of the process either with constraints [? ? ] 
or no constraints on the laser intensity [? ]. 

In this paper, we propose to do a step towards the ge- 
ometrical analysis of the control of dissipative quantum 
systems by beginning with the simplest system possible, a 
two-level system governed by the Lindblad equation. We 
determine control fields which minimize the total time of 
the process by applying the Pontryagin maximum prin- 
ciple (PMP) [????]. The maximum of the laser 
intensity is fixed to an arbitrary value. An increase of 
this value leads to a reduction of the duration of the con- 
trol. Note also that this cost functional seems particu- 
larly relevant in the context of a dissipative environment 
especially when the effect of dissipation is undesirable to 
reach the objective of the control. A second argument ex- 
plaining our choice of cost functional is of mathematical 
nature. The resolution of the PMP is particularly simple 
in this case and can be done analytically. The geometrical 
description of the time-optimal control is also well devel- 
oped especially on M. 2 (see [? ] for a recent overview) 
for affine systems with a drift term corresponding here 
to the dissipation. The main tools of this construction 
will be recalled throughout the paper. Finally, all this 
mathematical arsenal allows us to answer some physical 
questions such as the benefit that can be gained from 
dissipation. 

The paper is organized as follows. We first introduce 
the model for the Lindblad equation of a two-level sys- 
tem. Writing the density matrix in the coherence vector 
form [? ] and assuming some constraints on the control 
term, we restrict the problem to a control on a closed 
submanifold of R 2 with a true drift term i.e. an uncon- 
trollable term which cannot be eliminated by a unitary 
transformation. This drift term is due to the dissipation. 
We consider different cases with different initial and tar- 
get states that represent several situations of physical 
interest. We formulate in Sec. IIIII the PMP with a time 
minimum cost functional. In Sec. IIV1 we analyze the 
structure of the reachable set from each initial state and 
we point out the particular role of the fixed point of the 
free-dynamics in this structure. Section [V] deals with 
the construction of the optimal syntheses by solving the 
PMP. We recall that the PMP allows us to derive ex- 
tremal controls. Optimal trajectories are a subset of this 
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set which can be determined in a second step by a direct 
comparison of different extremal trajectories or by geo- 
metrical arguments. Conclusions and prospective views 
are given in Sec lVll Some technical calculations are re- 
ported in appendices lAl and IBl 



II. THE MODEL SYSTEM 

We consider a dissipative two-level quantum system 
whose dynamics is governed by the Lindblad equation. 
The system is described by a density operator p(t) which 
is a positive semi-definite hermitian operator acting on 
a two-dimensional Hilbert space Ti. The evolution equa- 
tion can be written as 

.dp 



at 



[H + uH x , p] + iC D (j>) , 



(1) 



where Ho is the field-free Hamiltonian of the system, Hi 
represents the interaction with the control field and Co 
the dissipative part of the equation. In the Lindblad 
equation, Cd can be written in the general form 

3 . 1 

(2) 



Since Ti[p] = 1, the density matrix p depends on three 
real parameters which can be given by the coordinates of 
the coherence vector [? ] : x\ = 25i[ ( 5i2], X2 = 23[/5i2] 
and X3 = P22 — Pii- From Eq. ([5]), one deduces that the 
coordinates Xi satisfy the following system of inhomoge- 
neous linear differential equations 

x\ = —Txi + u 2 x 3 

x 2 = -Tx 2 - uix 3 ,(7) 
£3 = (712 - 721) - (712 + 721)^3 + u x x 2 - U 2 Xi 

u\ and u 2 being two real functions such that u = m+iu 2 . 
As Tr[p 2 ] < I, we also have I < 1 which defines 

the Bloch ball. The dynamics is called either unital if 
712 = 721 i-e. the fixed point of the free dynamics is the 
center of the Bloch ball or affine otherwise [? ? ] . 

To simplify the study, we restrict the dynamics to a 
submanifold M C K 2 by assuming that the control field 
is real i.e. u 2 = 0. This hypothesis means that the control 
field is linearly polarized. With this choice of control, the 
last two equations of the system of Eqs. ([7]) are decoupled 
from the first one and the problem is reduced to R 2 . The 
system of Eqs. ([7]) then becomes 



E 

k=l 



[L k pL k \ - -{pL k \L k + L k \L k p)] 



where the operators L k are given by 



Li = 



■M = vr 



712, 721 and f are positive real constants describing the 
interaction with the environment. 712 and 721 correspond 
to population relaxations whereas T is the pure dephasing 
rate. We consider resonant fields, i.e. the frequency uj of 
the laser is equal to the energy difference between the two 
levels. In the RWA approximation, the time evolution of 
p(t) satisfies the following Redfield form of the Lindblad 
equation 



'at 



where u is the complex Rabi frequency of the laser field 
(the real and imaginary parts are the amplitudes of two 
orthogonal linearly polarized fields). T is the total de- 
phasing rate which can be written as 
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r = ^(712 



'721; 



(5) 



From Eqs. ([5]), we also notice that the requirement that 
Eq. (f5]) corresponds to a Lindblad equation implies the 
constraint T > 5(712 +721) or equivalently T > 0. Equa- 
tion ^ is written in units such that h — 1 . In the inter- 
action representation, Eq. ([4]) becomes 
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.(6) 



x 2 = —Tx 2 — 11x3 
x 3 = 7_ - 7 + x 3 + ux 2 

where 7_ = 712 —721 and 7 + = 712+721- The coordinate 
1 is set initially to 0. To simplify the notation, the index 



(8) 



has been omitted when confusion is unlikely to 
occur. Note that the analysis of the optimal control on 
M 3 is considerably more complex [? ] and goes beyond 
the scope of this paper. Equations © can be written in 
a more compact form 



x = F + uG , 



(9) 



with the vector x of coordinates (x 2 ,X3) and the two 
vector fields F and G defined by 



7- 



-ra;2 

- 7+Z3 



and G 



-X3 

%2 



(10) 



Pll\ r 

^kjal|y,'w straightforward calculation shows that the 
K fjbg|d/point of the free-dynamics is given for 7+ 7^ by 



x 2 = 

X3 = 1 - 

A 7+ 



(11) 



which is therefore either affine or unital according to the 
value of 7_. The irreversibility of the dissipation effects 
is reflected in the fact that the vector field F is a true 
drift term which cannot be eliminated by unitary trans- 
formations f? 1. 



III. PONTRYAGIN MAXIMUM PRINCIPLE 

We analyze the optimal control of this two-level system 
with the constraint of minimizing the total time of the 
control. We assume that the field u is bounded by 



\u\ < 1 



(12) 
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Equations ([8]) being linear, other bounds for u can be 
considered from a standard rescaling of the time and 
the dissipative constants T, j + and 7_. The Pontrya- 
gin maximum principle [????] is formulated from 
the following pseudo-Hamiltonian Hp 



H P = p- (F + uG) +p 



(13) 



where p = (p2,P3) € (M 2 )* is called the adjoint state and 
Po is a negative constant. We recall that the cost func- 
tion for the minimum time problem is equal to 1. This 
term is multiplied by the constant po m the Hamiltonian 
Hp. The Pontryagin maximum principle states that the 
extremal trajectories maximize Hp i.e. 



H max {yL, p) = Max\ u \<iH P (x.,v,u) 



(14) 



The coordinates of the extremal vector state x and of the 
corresponding adjoint state p fulfill the Hamilton equa- 
tions 



%2 
%3 



—Tx2 — ux 3 

7_ - 7+^3 + ux 2 



(15) 



x o 




and 



FIG. 1: Division of the manifold M by the sets Ca = A^^O) 
(in solid line) and Cb = A J3 1 (0) (in dashed and dot-dashed 
lines) for 7_ 7^ 0. The dashed and dot-dashed lines represent 
respectively anti-turnpike and turnpike arcs which are intro- 
duced in Sec. [V] The exterior circle in solid line corresponds 
to the limit of the Bloch ball in the plane (X2, £3). The small 
open circle indicates the position of the fixed-point of the free- 
dynamics. Numerical values are taken to be F = 1, 712 = 0.1 
and 721 = 0.3. 



(16) 



P2 = TP2 - Up 3 
P3 = 1+P3 + Up 2 



where u is here the extremal control given by Eq. (j 14(1 . 

We can now pass to the construction of the optimal 
syntheses. The construction begins with the introduction 
of two sets of points A7 1 (0) and A i3 1 (0) denoted Ca and 
Cp which divide M in different regions [? ]. A a and Ap 
are two functions from M to K defined as follows 



A A (x) = Det(F,G) 
A B ( x ) = Det(G, [F, G\) 



(17) 



where Det is the determinant of two vector fields and [., . 
their commutator. In our case, simple algebra leads to 



A^(x) = -Tx% + 7_£ 3 

A_b(x) = 2Tx 2 X3 — 2j + X 2 X3 + 7-£2 



(18) 



The major role of these two sets in the resolution of the 
optimal control problem will be detailed in Sec. [V] We 
can already say that Ca and Cp are responsible for qual- 
itative modifications of the optimal trajectories [? ? ]. A 
preliminary step thus consists in analyzing the structure 
of the sets Ca and Cp when T, 7+ and 7_ vary with the 
conditions 7+ > and T > From a formal point 
of view, this can be done by introducing the feed-back 
group [? ] but the dynamics being here bilinear [Eqs. 
flSJ], we can directly compute the sets Ca and Cp and 
their relative positions. 

If r 7^ 7+ then the set Cp corresponds to the union of 
the two lines 



and 



x 3 



-7- 



2r - 2 74 



(20) 



x 2 = , 



(19) 



with the restriction that \xs\ < 1. In the case V — 7+, Cp 
is only composed of the vertical line of equation x 2 = 0. 
The points (x 2 , x 3 ) of the Bloch ball that are solutions of 
the polynomial equation 7+x§ — 7^x3 — Tx\ — belong 
to Ca- In the case 7_ 7^ 0, Ca is therefore the union 
of two parabolas. This set is either above or below the 
line X3 — 2 r-2j + according to the signs of 7_ and of 
T — 7 + . For 7_ = 0, this set is reduced to the origin 
of the Bloch ball. Figure Q] displays the submanifold M 
together with the sets Ca and Cp for a given value of the 
parameters. We finally notice that the fixed-point of the 
free dynamics and the center of the Bloch ball belong to 
C A - 



IV. REACHABLE SETS AND 
CONTROLLABILITY 

We consider four different qualitative cases of control 
which allow one to study several physically relevant sit- 
uations : 

• Case (a) : Conversion of a pure state into a mixed 
state with a unital Lindbladian (7-/7+ = 0, T > 
7+ + 2). 

• Case (b) : Conversion of a pure state into a mixed 
state with a unital Lindbladian (7-/7+ = 0, 7+ — 
2<r< 7+ + 2). 
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• Case (c) : Purification of the completely random 
mixed state which corresponds to the center of the 
Bloch ball (7-/7+ = -1, T > 7+ + 2). 

• Case (d) : Conversion of a pure state into a mixed 
state with an affine Lindbladian (7-/7+ = —0.5, 

r>7+ + 2). 

The numerical values we have chosen for illustrations are 
given in Table HI Although the choice of the parameters 



TABLE I: Numerical values of the dissipative constants in 
arbitrary units. 





r 


712 


721 


(a) 


3 


0.3 


0.3 


(b) 


1.5 


0.3 


0.3 


(c) 


3 





1 


(d) 


3 


0.1 


0.3 



of Table |T] will become clearer in Sec. [V] some com- 
ments can already be made. This choice both depends 
on the structure of the sets Ca and C'b and on the char- 
acteristics of two particular extremals denoted X~ and 
Y — which start at the initial point, and correspond re- 
spectively to a constant control equal to -1 and 1. As 
detailed in the appendix [Bl the X— and Y— trajecto- 
ries are either pseudo-periodic or aperiodic according to 
the sign of the discriminant A = (T — 7+ ) 2 — 4 of the 
system of Eqs. ([8]). An exact resolution of the dynam- 
ics shows that the trajectory of the system is aperiodic 
if A > and pseudo-periodic otherwise. This point is 
summarized by the diagram [TO] of the appendix [Bj In 
Table HI we have chosen for three of the four examples T 
such that r > 7+ + 2 to simplify the local structure of 
the optimal synthesis around the fixed-point of the dy- 
namics. A pseudo-periodic trajectory is locally a spiral 
in the plane (2:2, £3) around this fixed point which makes 
the analysis more complex (see Sec. fV|) . 



A. Purity and limits of the dynamics 

Before analyzing the reachable set and the controllabil- 
ity of each example, we begin by some general comments 
about the dynamics. 

We first show that the field cannot locally compensate 
the effect of dissipation. The purity of the quantum state 



is defined by the function 2Tr[/> 2 ] — 1 = x\ + x\ 
states are thus on the unit circle of the {x2,x^] 
Simple algebra then leads to 



dt 



= -Tx 2 + 7.3:3 



. Pure 
plane. 



(21) 



points x where A^ (x) < then dTl jf J < and inversely 

if A A (x) > then > 0. The curve C A divides the 

plane (X2, X3) into a region where the purity of the state 
locally increases and a region where it locally decreases. 
On the boundary Ca, the purity is preserved. This point 
can be qualitatively understood as follows. We recall that 
the purity is equal to the square of the distance to the 
origin. The conservative vector field G is orthoradial (i.e. 
normal to radial vectors) for each point (£2,2:3) 7^ (0, 0) 
of the manifold. The dissipative vector field F does not 
modify the purity of the state if the radial component of 
F vanishes i.e. if F is parallel to G which is the definition 
of the curve Ca- 

We next analyze the fixed points of the dynamics when 
the field is on, which are defined by F + uG = 0. Since F 
and G are parallel, the fixed points belong to the curve 
Ca- The field-free limit point is the point of this line of 
maximum purity (see for instance Fig. [1]). This shows 
that the dissipation alone allows to reach the state of 
maximum purity. Inversely, one can ask if every point 
of the curve Ca corresponds to a limit point of the dy- 
namics. The answer is positive for a real non-bounded 
control u since the limits can be written 



— U7_ 

■ T 2 — r 7+ +« 2 

7_ 

^3 - 7+ +«vr 



B. Reachable sets 



(22) 



We begin by recalling some results of Refs. [? ? ] 
about controllability of dissipative systems. A quantum 
dissipative system of finite dimension whose dynamics is 
governed by the Lindblad equation is generically accessi- 
ble but not controllable. The accessibility property char- 
acterizes the fact that the system can be driven in every 
direction of the state space. Moreover, the concept of 
accessibility does not take into account the reversibility 
or the irreversibility of the process. The lack of control- 
lability is measured by the non small-time controllability 
of the system which illustrates the irreversibility of the 
dynamics. This kind of system is not small-time con- 
trollable because the field cannot locally compensate the 
effect of dissipation as shown by Eq. (f2"Tj) . The accessi- 
bility property can be checked by the computation of the 
dimension of the dynamical Lie algebra L of the system 
[? ? ] . We introduce for that purpose the density matrix 
p of coordinates (1,2:2, £3) which is given in the basis of 
the coherence vector and we rewrite Eqs. ([5]) in matrix 
form as follows 



p = CfP + uCcp , 
where Cf and Lq are 3x3 matrices given by 



(23) 



and we notice that this derivative does not depend on u, 
which completes the proof. A more general proof of this 
point is given in Refs. [? ? ]. We also point out that for 



C, 






^ 


f0 







-r 


and C G = 

















1 





.(24) 
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L is the Lie algebra generated by Lf and Lq. If 
r 7^ 7_|_, a direct computation shows that L is ci- 
ther isomorphic to gl(2) or to the semi-direct sum 
g[(2)(S)IR 2 of respective dimensions 4 and 6 for 7_ = 
or 7_ 7^ 0. The system is therefore accessible for f^7 + . 



We now determine the reachable sets from their respec- 
tive initial states of the four examples. The reachable sets 
can be constructed by the explicit construction of all the 
trajectories. We denote by 7?.(xo) the reachable set from 
xo. We first search for in this section the boundary of 
the reachable sets. Then in Sec. |V]we show that all the 
points inside the boundary are attainable. We consider 
for that the X~ and Y— trajectories starting from the 
initial point and from the field-free fixed point of the dy- 
namics. A qualitative change occurs when the X— or Y — 
trajectories cross Ca since the angle between the vector 
F(x) and G(x) changes its sign. The determination of all 
the extremal trajectories in Sec. IVl allows us to complete 
the construction. 

We apply these remarks to the cases (a) and (b) where 
the initial point has coordinates [xq. = 0, X3 — 1) and Ca 
is reduced to the origin (x 2 = 0, X3 = 0). The origin is 
attained asymptotically by the dynamics when t — ► +00. 
The Y '— and X— trajectories, which have the same initial 
and final points, are therefore global boundaries of the 
reachable set. This point is illustrated in Figs. [5^, and 
\Bp for the cases (a) and (b), where the reachable set TZ 
is in grey. 

For the case (c), the X— and Y— trajectories starting 
from (22 = 0,2:3 = 0) intersect asymptotically Ca- We 
also consider the X— and Y— trajectories starting from 
the field-free fixed point which define two new regions. 
As an infinite time is necessary to reach this point, these 
two trajectories do not belong to the reachable set. TZ is 
the union of all these regions as shown in Fig. We will 
check in Sec. E]by constructing the extremal trajectories 
that every point of this set is effectively attainable. 

For the case (d) the situation is more difficult to 
analyze as the X— and Y — trajectories starting from 
(x 2 = 0,Xs = 1) intersect transversally Ca- From this 
point of intersection, the two trajectories do not corre- 
spond anymore to the boundary of TZ. After the crossing 
of Ca, the X— or Y— trajectories originating from the 
preceding curves define two new regions of TZ. We finally 
consider the limit point of the field-free dynamics which 
can be attained as t — > +00 and the X— and Y — tra- 
jectories starting from this point. These two lines define 
two new boundaries which intersect asymptotically Ca 
but which do not belong to TZ. TZ is therefore composed 
of all the regions constructed above and is displayed in 
Fig. As before, only the determination of the ex- 

tremal trajectories will complete the construction. 
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FIG. 2: Reachable sets TZ (in grey) for the cases (a), (b), (c) 
and (d). The sets Ca and Cb are represented in dashed lines. 
The position of the field-free fixed point is indicated by the 
small open circle. In (c) and (d), the trajectories starting 
from the field-free fixed point do not belong to TZ (see text). 
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V. TIME-OPTIMAL CONTROL 

The application of the Pontryagin maximum principle 
is particularly simple at least locally in the case of time 
minimum control. The goal of this section is to solve this 
problem by constructing the optimal synthesis i.e. the set 
of all the optimal solutions starting form a given initial 
point Xo and reaching a point of TZ(xq). For time opti- 
mization, the optimal controls are composed of piecewise 
constant parts (u = ±1) and of singular controls. Some 
preliminary work in Sec. IV Al has to be done before de- 
termining the optimal solutions. 



A. Preliminary 

We define two functions $ and 9 [? ]. 
switching function, is given by 

= p.G = -P2X3 + P3X2 



$, called the 



(25) 



Using Eq. (|13j) . standard considerations of maximiza- 
tion of the pseudo-Hamiltonian Hp lead to the con- 
clusion that the extremal field at time t is given by 
u = sign($(t)) if <£>(£) ^ 0. A time to such that the con- 
trol changes sign (i.e. such that <{>(io) = 0) is a switching 
time. If $ vanishes on an interval [ioi^i] then the corre- 
sponding trajectory is called singular and is referred to 
as a Z— trajectory. An important point is the fact that 
the Z— trajectories lie in the set Cb and that the corre- 
sponding control u — <fi is singular i.e. different from 1 or 
-1.0 can be calculated by imposing that ^As(x(t)) = 
on [to, ti]. In our case, it can be shown that 



d . , .... . dA B 0A B 
5 A B (*(t))=0 = _«, + _ 



x 3 • 



One arrives after simple algebra to 

-x 2 j-(T-2 1+ )-2x2X 3 (7l-T 2 ) 
^ 2(r- 7+ )(z2-s§)-7-S3 



(26) 



(27) 



For the line X2 — of Cb, this leads to </>(x) = 0. For 



the line x 3 



2(r- 7 +) 



, we obtain 



(x) = 



7-( 7+ -2r) 
2(T- 1+ )x 2 



(28) 



The control is admissible if |</>(x)| < 1 which implies here 
the condition 



M > 



, 7-(7+-2r) 
1 2(r- 7+ ) 1 



(29) 



Not every Z— trajectory can be an optimal trajectory. 
We introduce to characterize this point the notion of 
turnpike and anti-turnpike curves [? ] . Let x ^ C'a U C b 
and the function /(x) = — sffej • A turnpike curve is an 
arc lying in Cb such that for every point x of this arc 



Aa(x) ^ and X(x) and Y(x) are not tangent to Cb- 
It is also assumed that X(x) and Y(x) point to oppo- 
site sides of Cb which define two regions tt x and il v . If 
/(x) > (resp. /(x) < 0) on fl y and /(x) < (resp. 
/(x) > 0) on Q x then the arc is a turnpike (resp. anti- 
turnpike) arc. The relation with the optimal synthesis 
can be stated as follows. Using for instance the clock 
form a (see appendix [XJ , it can be shown that the anti- 
turnpike trajectories are not optimal. Figure [1] displays 
the turnpike and the anti-turnpike curves for particular 
values of the parameters. 

The switching function 4> is a powerful tool that gives all 
the informations on the extremal trajectory. However, a 
global analysis of the extremals using only $ is difficult 
because all the initial values (p20jP3o) have to be tested. 
A more global point of view is given by the function 9 
which has the advantage not to depend on p. 

Let v be the vector of coordinates (2:2, x'3). Deriving 
Eqs. © with respect to time, it can be shown that v 
satisfies the following system of equations 



V2 
V3 



~Tv 2 - uv 3 
-7+V3 + uv 2 



(30) 



We next introduce the vector v(t). At time t, let us as- 
sume that the dynamics reaches the point (xz(t) , xs(t)) 
and that the vector field G at this point is given by 
G(x2(t),X3(t)) = G(t). v(t) is defined as the solution 
at time of Eqs. (|3"0")) that at time t is equal to G(t). 
The dynamics is thus propagated backwards during the 
time t. 9 is then defined as 



0(t)=arg(t;(O),t;(t)) , 



(31) 



where the angle is measured counterclockwise. The def- 
inition of 9 originates from the fact that the product 
p(t) • v(t) is a constant [? ]. We thus have 



*(t) = p{t) ■ G(t) = p(0) ■ v(t) 



(32) 



The function 9 allows to determine some properties of the 
extremal trajectories. We briefly recall these points here. 
The reader is referred to [? ] for rigorous definitions and 
complete proofs. It can first be shown that 



sign(9) = sign(As) , 



(33) 



which means that the zeros of 9 are located on Cb- A 
switch can occur if 9 > and 9 > or 9 < and 9 < 0. 
In addition, the variation of 9 between two switches or 
between a switch and a singular control is equal to 
modulo 7r. This latter point can be easily understood 
from Eq. (f3"2"| . Other properties of 9 will be detailed 
during the construction of the optimal synthesis. 



B. Time-optimal syntheses 



The parts lV B 11 IVB21IV B 3l and lV B 4l are rather tech- 
nical and describe the way to obtain the optimal synthe- 
ses. Section IV B 51 details conclusions on the role of the 
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dissipation that can be gained from the resolution of the 
time optimal control. 

1. Case (a) 

In the case (a), the optimal trajectories are either bang 
or bang-bang. A bang trajectory is a trajectory associ- 
ated to a single value of the control u = +1 or u = — 1. 
A bang-bang trajectory is the concatenation of an X— 
and an Y — trajectories. We denote by X * Y such a con- 
catenation where the Y— trajectory comes first. Here 
the maximum number of switches is thus equal to 1. The 
elimination of extremal trajectories selected by the Pon- 
tryagin maximum principle has been done through the 
clock form a (see appendix [Al for details) and the sym- 
metry of the diagram with respect to the line X2 = 0. We 
recall that the clock form can only be used for trajecto- 
ries belonging to one of the four quadrants defined by Cb 
and which do not cross Ca- Figure [3] displays the optimal 
synthesis for this problem and the evolution of the angle 
along the Y— trajectory starting from the initial point 
(0, 1). The plot of 8 shows that the switches are always 
permitted along the X— or the Y— trajectory. Note also 
that the form of the curve representing 8 is related to the 
sign of A. More precisely, if T > -f + + 2 then 8 is a mono- 
tonically decreasing function and if T < j + — 2 (case not 
treated here), 8 starts increasing, passes through a maxi- 
mum when the trajectory crosses Cb and then decreases. 
The intermediate case 7+ — 2 < T < 7+ + 2 corresponds 
to the case (b). Using the clock form, it can be shown 
that only one switch is possible for a trajectory in a given 
quadrant. This point can also be determined from more 
general considerations detailed in Ref. [? ]. The rest of 
the optimal synthesis is deduced from the symmetry of 
the problem. The line xi = is called an overlap curve 
denoted by K as it is the locus of points reached by two 
optimal trajectories. We finally notice that the dissipa- 
tion alone is not used here by the control and therefore 
cannot help accelerating the control. We recall that the 
singular control on the vertical line of Cb is given by 
u = 0. 



2. Case (b) 

The situation is a little more complex in the case (b). 
For X3 > 0, the synthesis is similar to the one of case (a) 
i.e. the trajectories are either bang or bang-bang. The 
function 8 is here periodic, the times where 8 vanishes 
correspond to the crossing of Cb- This function tells also 
us that the initial Y — and X— trajectories cannot switch 
in their respective second quadrant (i.e. the second quad- 
rant they go through) since (t) > and 8{t) < 0. By 
symmetry with respect to the line 22 = 0, one deduces 
that theses curves are optimal up to this line. A singular 
control along the horizontal line of Cb is optimal from 
the point of intersection between the initial X~ or Y — 



Y 






X 















-0.4 -0.2 0.2 0.4 




t 



FIG. 3: The top figure represents the optimal synthesis for 
the case (a) (see text). The dashed lines indicate the locus 
Cb and the small open circle the fixed point of the dynamics. 
The bottom figure displays the evolution of 9 as a function of 
the time for the Y— trajectory starting from the point (0, 1). 



trajectories and this line. Since 7_ = 0, the singular 
control is given by <fi = 0. From this singular line origi- 
nate optimal trajectories with control u = ±1. Using the 
clock form a, it can be shown that these trajectories can 
not switch again. To summarize, in this second part we 
have constructed optimal controls of the form Y * Z * X, 
Y*Z*Y, X*Z*X and X * Z * Y . This optimal synthesis 
is represented in Fig. 3J Here, one sees that the dissipa- 
tion alone (with u = 0) accelerates the control to reach a 
point near the origin only along the horizontal direction 
and for X3 = 0. We finally give in Fig. [5] an example of 
the comparison of two extremal trajectories. The time 
to reach respectively the points A and B from the initial 
point (0,1) by the Y — and X— trajectories is the same. 
We would like to attain the point C from A or B. We use 
either the A— trajectory from B or a concatenation of a 
Y — and a Z— ones from A. The clock form a cannot be 
used here since the trajectories belong to two different 
quadrants. We therefore consider the symmetric image 
AD of the curve BD with respect to X2 = 0. The two 
extremals are thus in the same quadrant and a can be 
used. We conclude that the curve Y * Z * Y is optimal 
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FIG. 4: Same as Fig. [3] but for the case (b). The clotted- 
dashed line represents the singular trajectory S. The small 
insert is a zoom of the optimal synthesis near the origin. 
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X 2 

FIG. 6: Same as Fig. [4] but for the case (c). 



3. Case (c) 

In this case, the singular line X2 = is optimal. The 
singular control is equal to zero. Switches can occur from 
the initial X— and Y~ trajectories but they do not lead 
to optimal trajectories. Inversely, X— and Y — trajecto- 
ries originating from S are found to be optimal. When 
two extremals cross Ca, ol cannot be used and a direct 
numerical comparison is then performed. In this case, 
the function 9 tells us that these curves do not switch. 
Trajectories originating from S that have a switch are 
therefore not extremal. We finally point out that the 
dissipation is beneficial for the control and decreases the 
duration needed to purify the system. The optimal con- 
trols are here of the form X—, Y-, X*Z or Y*Z. Figure 
|5] displays the optimal synthesis for this problem. 




FIG. 5: Use of symmetry to determine the optimal trajectory 
(see text). 



to reach the point C. 



4- Case (d) 

The case (d) is the more complex one and corresponds 
roughly to the composition of cases (b) and (c). The 
difficulty lies in the global structure of the control or, 
in other words, in gluing the two preceding local analy- 
sis. For X3 > 0, the optimal synthesis is similar to the 
cases (a) or (b) with bang or bang-bang trajectories. The 
bottom of the optimal synthesis from the point of inter- 
section of the initial X~ and Y — trajectories is similar 
to the case (c). 

We now describe the central part of the synthesis. The 
horizontal singular line of Cb does not correspond to 
a singular trajectory since |<^(x)| > 1 which is a non- 
admissible control. Consider the first point of intersec- 
tion of Cb with the X— or Y — trajectory. Following [? 
], we know that a switch curve C originates from this 
point. To determine the exact locus of C, we use the 9 
function since this function varies by modulo it between 
two switches. C is constructed numerically. By defini- 
tion of C, every trajectory which crosses C switches on 
C. We have found that C, Cb and Ca intersect in one 
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X 2 

FIG. 7: Same as Fig. [4] but for the case (d). The dotted line 
represents the switch curve C. 



point, the origin O. 

Since the line x 2 = is turnpike for 2=. < ^ 3 < 0, we 
can ask if this singular trajectory is optimal i.e. if we 
can have a local optimal synthesis of the form given by 
Fig. [5] To answer this question, we use the switching 
function For x(t) e CUS, $(t) = i.e. the vectors 
p(t) and G(x(i)) are orthogonal. Since the direction of G 
is known (G is orthoradial), one can deduce the direction 
of p(t). Let zi and z 2 be two points belonging respec- 
tively to G and 5*. The vectors G(z) associated to these 
points are schematically represented in Fig. [51 We let 
now the states Zi and z 2 go to (0, 0) and we determine 
the directions of the different adjoint states. We recall 
that the Pontryagin maximum principle states that p is 
a continuous function which does not vanish. When Zi 
goes to (0,0), one deduces by a continuity argument that 
is vertical in O. When z 2 goes to (0,0), the limit di- 
rection of p 2 is given by the switch curve G. To respect 
the continuity of p, one sees that G has to be tangent to 
the line x 2 — in O. Due to the complexity of analytical 
calculations, we have checked numerically that this is not 
the case. The singular line for x% < is therefore not op- 
timal. The trajectories of the form Y*X*Y or X * Y * X 
are thus optimal up to x 2 = 0. In addition, when the ini- 
tial X— and Y— trajectories cross Ca, the angle between 
the vectors -F(x) and G(x) changes its sign. New optimal 
trajectories originate from this point of intersection and 
correspond to two new regions of the reachable set. 



5. Qualitative conclusions on the dynamics 

From the results obtained in the preceding sections, 
some qualitative conclusions can be made with respect 
to the dissipation effect on the time optimal control of 
the dynamics. The dissipation is not undesirable when 
the dissipation allows to purify the system [cases (c) and 
(d)] and help accelerating the control. In contrast for the 
conversion of a pure state into a mixed state, the dissi- 




FIG. 8: Possible optimal synthesis around the origin. 



pation alone increases the duration of the control and its 
effect is not beneficial for the control [cases (a) and (b)]. 
The example (d) summarizes well the situation. As long 
as the purity of the state decreases, it is advantageous to 
use a control field but when the purity starts increasing 
the dissipation alone becomes more efficient. 

Finally, we point out that a visual inspection of the op- 
timal syntheses allows to check the form of the reachable 
sets given in Sec. IIVBI 



VI. CONCLUSION 

We have investigated the time-optimal control of a 
two-level dissipative system. For different initial and fi- 
nal points, we have constructed the optimal syntheses 
which allow us to rigorously conclude on the beneficial 
role of the dissipation in order to reduce the duration 
of the process. An open question is the generalization 
of the present study to more complex quantum systems 
with for instance two control parameters or three lev- 
els. The determination of the optimal syntheses for these 
systems requires more sophisticated mathematical tools 
than those used in this paper [? ]. In addition, such 
an analysis could also be performed on other problems of 
quantum control. We point out that the question of mea- 
surement in the context of geometrical control theory is 
particularly interesting. Measurements can be viewed as 
a non-unitary control which allows to enlarge the reach- 
able sets. For instance, we can consider systems not com- 
pletely controllable by unitary control but which become 
controllable if measurements are used. A simple example 
is given by a finite number of levels of the harmonic os- 
cillator with a constant dipolar interaction only between 
adjacent energy levels [? ]. From a formal point of view, 
this problem raises the challenge of the precise formula- 
tion of the Pontryagin maximum principle in this context 
which is to our knowledge an open question. We notice 
that all these geometrical tools would allow us to an- 
swer some physical relevant questions such as the time 
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at which measurements have to be done to minimize the 
duration of the control. 



APPENDIX B: ANALYTICAL DETERMINATION 
OF THE DYNAMICS 



APPENDIX A: THE CLOCK FORM 

We derive in this section the expression of the clock 
form denoted a [? ]. By definition, the clock form is a 
1-form which fulfills the following conditions 



(Al) 



a(F) = 1 
a(G) = • 

A solution of this system exists except on the set Ca 
where F and G are collinear. If we write a as a = 
a 2 dx 2 + a 3 dx 3 then simple algebra shows that a 2 and 
a 3 are solutions of the system 



a 2 (-rx 2 ) + a 3 (7- 
a 2 x 3 = a 3 x 2 



7+^3) = 1 



We obtain that 



—X3 



(A2) 



(A3) 



rx|-7_x 3 +7 + 3:| 

From the 1-form a, we can define the 2-form da which is 
given by 



da 3 O0.2 

da = (- — )dx 2 A dx 3 

0x2 0x3 



and reads after some calculations 



da — 



2Tx 2 x 3 + 7_x 2 - 27+2:2X3 



[Txt 



l-x 3 +j+xl} 2 



dx2 A dx 3 



(A4) 



(A5) 



If we write da as da — g{x 2 ,x 3 )dx2 A dx 3 then one sees 
that g{x 2l x 3 ) — on Cb and that the function g has 
a constant sign in the regions delimited by the lines of 
Cb- This point is displayed in Fig. [9] As suggested by 
its name, the clock form is a form which allows one to 
determine the time taken to travel a path and to compare 
the extremals. Let 7 be a path in the plane (x 2 , x 3 ) which 
does not cross Ca and T the time of travel along 7. We 
have 

/ a = [ a(x)dt = [ a{F)dt = T . (A6) 
J j Jo Jo 

We consider now two paths 71 and 72 starting and end- 
ing at the same points and respectively associated to the 
durations T\ and T 2 . One shows by using the Stokes 
theorem that 

T x - T 2 = / a- a= da , (A7) 

J71 ^72 •* D 

where D is the surface delimited by 71U— 72. For paths 71 
and 72 which lie in one of the four quadrants defined by 
Cs, we can straightforwardly determine the time-optimal 
trajectory. 

*X 



g<o 



g>o 



In this section, we solve analytically the dynamics of 
the system given by the following system of differential 
equations 



a?a = — Tx 2 — ux 3 

x 3 — 7_ - 7 + x 3 + ux 2 



(Bl) 



We assume that u = ±1 and we denote by u = e the con- 
trol term. Combining Eqs. (|B1[) . one arrives to equations 
that depend on only one variable X2 or x 3 



±2 + (r + 7+)a; 2 + (1 + r 7+ )x 2 + £7- = 

± 3 + (r + 7+)x 3 + (1 + r 7+ )x 3 - r 7 _ = 



(B2) 



Equations (|B2j) are linear inhomogeneous differential 
equations of second order whose discriminant A is given 
by 



a = (r - 7+ ) 2 - 4 



(B3) 



A qualitative change of the solutions is expected accord- 
ing to the sign of A, i.e. a transition from a pseudo- 
periodic solution for A < to an aperiodic solution for 
A > 0. This point is summarized in Fig. [TO] Note that 
the parameter 7_ plays no role in the structure of these 
trajectories but modifies their limits points. Simple al- 
gebra allows one to calculate the exact solutions. For 
instance for A > 0, we obtain 



f x 2 (t) = e^ r +^ 2 [x2ocosh(lt) - (x 2 oi - x ^ + f x 30 )sinh0 
\ X3 (t) = e-( r +^)*/ 2 [a; 3 oCOsh(|i) + (^^x 30 + ^x 2 o + ^)sin 



where (x 2 o,x 3 q) is the initial point of the dynamics and 
S = \/~K. 
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FIG. 10: Sign of the discriminant A as a function of the 
parameters T and 7+. The parameters are in arbitrary units. 
The zone in grey is excluded for Lindblad dynamics. 



